Absence of Quantum Metallic Behavior in Disordered Granular Superconductors 
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We examine the idea, postulated by Phillips et al., that a finite resistivity in T ^ limit in 
disordered granular superconducting (SC) films is explained as a consequence of the absence of 
phase stiffness in a phase glass (PG) peculiar to granular systems. It is found that, in spite of the 
absence of static phase stiffness, a coupling between the nonzero PG order and the ordinary SC 
fluctuation makes the conductivity divergent. However, an actual drop of resistance is argued to 
occur due to another SC glass ordering, induced by the precursory PG fluctuation, corresponding 
to the vortex-glass transition in a nonzero magnetic field. 
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Study of the phase diagram in disordered granular su- 
perconductors at low temperature (T) is important in re- 
lation to understanding the superconductor-insulator (S- 
I) transition behaviors. Recently, Dalidovich and Phillips 
have argued that, in a two dimensional (2D) phase 
glass (PG) expected to be induced in low T limit by a dis- 
ordered granular structure, the resistance is finite. Ob- 
servations of a quantum metallic behavior 0,0,0] do not 
seem to be directly explained within available theories 

Qi m for homogeneously disordered superconducting 
(SC) materials. Further, it is found that their analysis 
can be trivially extended to the 3D case and the case with 
nonzero magnetic field {H ^ 0) and leads to a similar 
quantum metallic behavior. Since resistive data showing 
a metallic resistance curves flattening upon cooling are 
available in quasi 2D systems in H ^ 0, their pro- 
posal, if correct, might become a correct description of 
the corresponding phenomena. On the other hand, it has 
been argued that, in contrast to the homogeneously- 
disordered case 0, the resistivity drop at low enough 
(but finite) temperatures in disordered superconductors 
consisting of clustered SC islands (i.e., SC grains) should 
occur even in a wide field range above an averaged (mean 
field) upper critical field -ffc2(0) . Therefore, effects of 
disordered granular structure on quantum resistive be- 
haviors are not theoretically understood well at present. 

To examine this problem, we use essentially the same 
model as in Ref.^H expressing a random Josephson junc- 
tion array with an on-site charging energy. 
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where the pair < i, j > denotes a nerest-neighbor pair, 
and 



So = 
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Jvir) = j;Ar) - J^, exp [i(A^f + <5Ay(r))], 

S±,{t) = exp(±i0,(r)). 
The gauge field consists of the static component 



for an applied uniform magnetic field H and the dynam- 
ical fiuctuation (5A(t) introduced for deriving the linear 
responses. The quenched disorder in the system is in- 
corporated in a randomness of with nonzero mean Jq 
{J'ij — Jo > 0) for any < i, j > and with a Gaussian 
distribution. Only for a formal justification of the mean- 
field approach on the PG ordering, the Gaussian distribu- 
tion on Jij — Jq will be replaced by an infinite-ranged one. 
Then, after replicating the action and introducing PG 
and SC order parameter fields, <7°''(ti, r2) — (g'"'(r2, ti))* 
and ■i/'f (''■) , the random-averaged free energy is given by 
lim„^+o(-^ — !)/"■, where 
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Zq is the replicated partition function for Sq, 

^os{4>.q)^^j dnj dT2 ^|g('^'')(ri,r2)|^ 
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and the bracket ( )o denotes the ensemble average by the 
charging energy part So- In eq.(4), a and b are replica 
indices, d is the spatial dimension, • D* denotes the 
Laplacian on the gauge- invariant gradients defined on the 
cubic lattice (for d = 3), and the factor N (the sys- 
tem size) in the first and second terms arises due to the 
replacement into the infinite-ranged model. The gauge 
fiuctuation SA — 6 Ax applied externally is assumed to 
be small, and its site dependence was neglected because 
spatially uniform linear responses are considered below. 
Performing the cumulant expansions with respect to g'"'') 
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and is obtained in a form of a Landau free en- 

ergy functional. Further, as in Ref.fTll H^ . rewriting 
g("''Hn,T2) as Q^''^\ti,t2) - CSa,bS{Ti - T2) in order to 
delete the term / dri J dT2|Q^°''^(Ti, T2)P, we finally ob- 
tain the following effective action 



zero-dimensional case) with Cg — > will be assumed be- 
cause extending to higher dimensional and H ^ cases 
can be trivially performed. Then, the variational equa- 
tion — hm„^+o"~^5Z"/9(5(°'')(Ti, is reexpressed 
by the following three equations: 
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K , 1(^2 + r) - kD^ + u{q + P'^ ^-') " w^Gd{uj) = 0, 
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For a moment, a |^/;|^ term with positive coefficient will 
be neglected. The coefficients in eqs.(5) and (6) can be 
calculated from correlation functions on S±{t) on the 
basis of the local action Sq, and, except r and r^, all of 
them are always positive . For instance, the coefhcient 
is postive as far as Jq > 0. The above expression of 
the action is of the same form as in Ref . jT^ . 

Now, let us examine the mean field solution on Q"'' 
by taking account of the presence of SC (-0-) fluctua- 
tion. In iJ 7^ case, this treatment may be appro- 
priate above -ffc2(0) and in low enough temperatures 
[iJl- Note that there is no reason why the SC fluc- 
tuation is negligible in the mean field analysis of the 
PG order because the tp-Re\d in eq.(7) couples to the 
Q- field in a bilinear form (t/i^"))*?/!^^). As shown be- 
low, including SC fluctuations is essential to obtaining 
a correct result of conductivity. Following Ref. [TIL , 
the mean field ansatz Q^"''^(ti, T2) = q{l — Sa,b) + (9 + 
T,^^oDuje^''^'-''^^''^^)Sa,b will be used together with 

4"'\t) = /3"^Ec^V'i''H'^)e"'"^- The assumption of 
replica symmetry in w = terms is sufficient for the 
present purpose of showing a divergent conductivity be- 
cause a breaking of replica symmetry would be accompa- 
nied by an independent parameter such as a coefhcient 
of a quartic term on Q"^. Further, only for the con- 
venience of presentation, the on-site T/j-fluctuation (i.e.. 



-«;^Gd(0) = 0, 
(9) 

(10) 



2Kf3^q(7! - q) + w^,GodiO) ^ 0. 

Here, the SC fluctuation propagator G'(°'')(cj) = /J^^ < 

> is given by G^'^^'^uj) = 5a,b{^ " 
5u.fi)Gd{uj) + (5„,o('^a,fcGrf(0) + (1 - Sa,b)God{0)), where 
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and Ag = j3{q — q). Noting that, when g > 0, eq.(lO) 
becomes 



Ag=(^) |:(d^(Ag))^ 



(12) 



we easily find that the only physically meaningful solu- 
tion of the PG order parameter is given together with 
eq.(12) by 



= -Aq - K 



u ^ {K{Aq)'^ + w^pd^ — K -^r), (13) 
1 + Kc^w^d'^/t 



tK-^d^^{w^/t)3 



The above form of D^^ is valid up to 0(|ct;|). In —f 
limit where -0 and Q fields are decoupled, the above 
solution reduces to the pure mean-field solution [3, 0, 
IT^ with q > and -0 = 0. The crucial point is that 
Du)^o ^ in the presence of the ■(/'-fiuctuation. Situation 
is similar to the Ising spin-glass case in a magnetic field 
[Tl|. Actually, it will be recognized that, in eqs.(8) to 
(10), the fluctuation propagator G'"''-' plays similar roles 
to an external magnetic field in the spin-glass problem. 
Further, it will be clear that eq.(13) is also valid in higher 
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dimensional case and H case if the expressions of d^, 
where m > is an integer, are appropriately replaced. 

Let us turn to examining the conductivity a in terms 
of Kubo formula. Contributions to a arise from J^a and 
the \Di^pi\'^ term in Due to the |a;|-dependence in 
D^j, the latter contribution leads, as in Ref.P|, to a finite 
contribution 1] to in T — > limit. The expression of 
J-'a clearly implies the absence of static phase rigidity. 
However, it does not imply a finite conductivity because 
oi pq ^ D^^Q 7^ 0. The PG contribution to cr arising 
from Ta, 



crpG(iw) 
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dT, Sin - T3)|Q("'')(n,T2)P - IQ^'^'Ht, - T3)|^ 



is essetially the same as eq.(14) in Ref. j2i and corresponds 
to the sum of a^^^ and cr^^^ defined in Ref.0- The first 
term cr'^^ results only from the static components q and q 
and thus, should vanish consistently with the absence of 
static phase rigidity. Actually, it was verified 13] through 
a delicate manipulation of an analytic continuation. We 
will focus rather on cr*^^^ , which is given by 
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We are not interested in the first two terms under the wi- 
summation which will lead to a finite contribution to a^^^ . 
The last term in the bracket, —2qDo, is positive at low 
T (= /3~^) limit and when q > because q :^ q— 0(T). 
Since, according to eq.(13), the second two terms in the 
bracket also give only a 0(T) correction to the above 
term, a divergent conductivity in low T limit results from 
eq.(15), as in the Meissner phase, as far as the PG order 
is present (i.e., q > Q) as a, consequence of a nonzero 
Jo (X . 

Finally, we briefly show that this divergence of con- 
ductivity due to the PG order is preceded by that due to 
another SC glass ordering induced by the fluctuation of 
PG order parameter Q"''. This SC glass corresponds to 
the vortex-glass 0| studied so far for nongranular mate- 
rials in 7J 7^ 0, and we closely follow below the treatment 
in Ref. fisj . Let us relax the infinite-range approximation 
for the disorder part (cx J^) of the action to incorporate 
a gradient term (VQ"'')^. Using results in Ref. 11] in the 
q — Q case and integrating over the Q^^-fluctuation in 
the Gaussian approximation, one finds that the nonlin- 
ear terms in ^l) of take the form jlSj 
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where -u^ > 0, and 0| 
C(k, wi,a;3, W2) = 



-G(k, ivi , a;3)G(k, , ^^4) 



1 -I- M/3-iE'^(^'^ + '^i'^ + ^3) 
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where 



Pu)\ .U12 (k) 
normalized. The 



is the Fourier transform of 
and the length scale was properly 
-C^ term in eq.(16), due to its ten- 
dency enhancing with decreasing f, dominates over the 
original term and suggests a great enhancement 9] of 
quantum SC fluctuation due to the granular structure. 
The onset of PG order is signaled, within the Gaussian 
fluctuation, by f +0, while the SC glass ordering is 
signaled by a divergence of the glass susceptibility 
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If 



expressing Xsg by ladder diagrams like Xa 
its irreducible vertex / 
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'Y^j^^i its irreducible vertex / is proportional to 
the vortex-pinning energy term, i.e., the minus second 
term of eq.(16), with vertex corrections due to the ij}- 
interaction term (the first term of eq.(16)). Let us re- 
strict ourselves to the lowest Landau level (LLL) of the 
'(/'-fluctuation and assume the presence of PG order in low 
T limit (i.e., r(T -> 0) ^ -f-O). Since the V- fluctuation 
in LLL is noncritical at low T limit, a f-dependence of / 
primarily determines Xsg- If the interaction vertex cor- 
rections are neglected, we find that / in LLL is divergent 
upon cooling like f"^/'', suggesting that the SC glass or- 
dering occurring prior to the PG ordering induces a di- 
vergence of conductivity. Hence, it is important to verify 
whether this divergence is suppressed or not by the inter- 
action vertex corrections. If calculating this vertex cor- 
rection in a consistent way with the Hartree approxima- 
tion [TsIITgI I. the strength (x r~^/* of the resulting renor- 
malized vortex-pinning energy is modified in the manner 



VqyG{qyy,Q,Q)l3 ^Y.^G{qyy,uj,uj) 



(1 



(19) 

_ 2 

where Vq^ = e ^y' , and qy is an external wave num- 
ber carried by the i/i-fluctuation propagator. Further, 
the factor in the bracket of eq.(17) was neglected be- 
cause it plays no roles in f — > limit, and the posi- 
ive coefficients w^, u, and t were set to be unity be- 
:ause their detailed values are unimportant here. Fo- 
cusing on the T ^ limit, one easily verifies that the 
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vertex correction merely brings an enhancement factor 
~ (ln[min(g^, Vf)])^. Thus, the result obtained without 
the interaction vertex correction is essentially unaffected, 
and, even in H > _ffc2(0), the SC transition (vanishing of 
resistivity) in nonzero fields is expected to be driven by 
the vortex-glass ordering induced by an effective vortex- 
pinning effect. 

This result, independent of the dimensionality of sys- 
tems, is apparently similar to but different from the main 



conclusion in Ref. ^3] ■ In Ref. ^3 1 an upper limit of the 
SC transition field at T = is obtained above i?c2(0) 
within the range where a random XY model similar to 
eq.(l) is applicable. On the other hand, the critical PG 
fluctuation inducing the SC ordering in the present model 
exists in any H as far as the granular model is valid and 
the SC fluctuation is present. Namely, an upper limit of 
T = SC transition field in the present case is roughly 
given by the field at which the model eq.(f) breaks down. 
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